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Abstract 

We analyze the quasinormal modes of the D2-D8' model of 2+1-dimensional, 
strongly-coupled, charged fermions in a background magnetic field and at non- 
zero density. The model is known to include a quantum Hall phase with inte- 
ger filling fraction. As expected, we find a hydrodynamical diffusion mode at 
small momentum and the nonzero-temperature holographic zero sound, which 
becomes massive above a critical magnetic field. We confirm the previously- 
known thermodynamic instability. In addition, we discover an instability at 
low temperature, large mass, and in a charge density and magnetic field range 
near the quantum Hall phase to an inhomogeneous striped phase. 
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1 Introduction 

Understanding the behavior of strongly interacting electrons in condensed matter 
systems is a long-standing and difficult problem. Among gapless systems with unbro- 
ken translation and charge conservation symmetry, Landau-Fermi liquids are the one 
well-understood class; much less is known about compressible systems that are not 
continuously connected to free fermions. Non-Fermi liquids, such as strange metals, 
which have Fermi surfaces but exhibit very different low-energy behavior, have been 
a hot topic recently, but they remain essentially mysterious. There is also the con- 
troversial possibility of compressible states with unbroken symmetry that are neither 
Fermi liquids nor non-Fermi liquids because they lack Fermi surfaces altogether. 

Studying the possible compressible states of fermionic matter would appear to 
be a problem amenable to gauge/gravity duality. Holography provides tractable de- 
scriptions of strongly-coupled systems and elegantly avoids the fermion sign problem 
typically plaguing calculations at nonzero density. 

We will take a top-down approach, with the goal of constructing a broad class of 
related string theory models and investigating both their common, universal features 
and the range of possible phenomena. Perhaps the most tractable and well-studied 
such group is the family of Dp-Dq brane intersection models, which includes, for ex- 
ample, the supersymmetric D3-D7 jlj, the related non-supersymmetric D3-D7' [2], 
the Sakai-Sugimoto [3], and the D3-D5 models 13], as well as many other variants 
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in 1+1, 2+1, and 3+1 dimensions]^] These brane systems describe an unusual com- 
pressible phase, the holographic quantum liquid; despite compelling arguments that 
a Fermi surface ought to exist [6] , various attempts^] have failed to find convincing 
evidence of one [7-11 . It is currently an open problem whether the Fermi surface 
is really absent or whether it is fractionalized or otherwise hidden from the view of 
gauge invariant observables. 

One interesting new member of this family of holographic quantum liquids is 
the D2-D8' system jTi], a holographic model for 2+1-dimensional, strongly-coupled 
fermions. This system, a close cousin of the D3-D7' modej^] |2|, consists of a stabilized 
probe D8-brane in the background of a stack of D2-branes such that #ND = 6. 
Unlike the D3-D7' model, the D2-D8' system exhibits hyperscaling violation, with an 
exponent 6 = — | 
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This is simply a consequence of the fact that the coupling 
constant is running trivially Indeed, one can trace back the conformal behavior in 
effectively higher dimensions by switching to a dual frame 18 . At nonzero charge 



density and magnetic field, the spectrum contains a gapped quantum Hall state with 
filling fraction one, whose excitations were analyzed in 19 and where, in particular, 
the existence of the magneto-roton as the lowest excitation mode was established. In 
this paper, we focus instead on excitations of the gapless, metallic states. 

We find a pattern of low-lying quasi-normal modes qualitatively similar to that 
found in other Dp-Dq systems [8 ^lT|[20]j25] . At long wavelength the longest lived 
mode is, as expected, a diffusion mode associated with charge conservation. At shorter 
wavelength, crossing over from the hydrodynamic to collisionless regimes, the diffusion 
mode merges with another purely dissipative mode, and the holographic zero sound 
mode appears. Turning on a magnetic field rearranges the spectrum of quasinormal 
modes and, when sufficiently large, gives a mass to the zero sound. 

Unlike the stable quantum Hall state, the metallic state is unstable in large regions 
of parameter space. We find two types of instability: a thermodynamic instability, 
which occurs next to the quantum Hall state, and spatially modulated instability at 
low temperature and for magnetic fields roughly of the same size as the charge density. 
The first of these was observed in 14 as a negative charge susceptibility. Here, we see 
this same instability manifest in the linear fluctuation analysis as a diffusion mode 
with negative diffusion constant. 

The modulated instability is induced by a nontrivial term in the D8-brane Chern- 
Simons action and is in the family of instabilities seen first in 26 . Instabilities at 



nonzero momentum due to Chern-Simons terms turn out to be quite generic and have 
been found recently in a variety of holographic systems 22 23 27 . We have not yet 



1 For a partial list, see JHj. 

2 Bottom-up investigations of holographic matter have uncovered many phenomena which are 
intimately related to the existence of the Fermi surface. Two recent examples are Friedel oscillations 
and Landau levels 
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3 Also, see 
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16 



13 



for earlier developments leading up to these types of models. 
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found the inhomogeneous ground state which follows from this instability, because 
the worldvolume fields on the D8-brane are fully coupled. We anyhow expect that it 
will be a striped phase with all quantities, such as spin, charge density, and transverse 
currents, spatially modulated. 

Although our results are in many ways qualitatively similar to those found in 
the D3-D7 and D3-D7' models, there are some differences worth highlighting. The 
D3-D7 model is supersymmetric and free of instabilities [7]. The D3-D7' system, on 
the other hand, suffers an analogous modulated instability to the one found here. In 
D3-D7' model, increasing the charge density exacerbated the modulated instability, 
while adding a magnetic field or increasing the mass improved stabity, in line with 



expectations from the original supergravity model 26 . Surprisingly, we find that 



the D2-D8' system is unstable near the quantum Hall state, that is, when the charge 
and magnetic field are roughly equal. When either the charge or magnetic field is 
much larger than the other, the homogeneous state becomes stable. Also, in contrast 
with the D3-D7' result, increasing the mass enhances the instability. Furthermore, 
the thermodynamic instability observed earlier in 14 near the quantum Hall state 
is present only in D2-D8' model and is found to be continuously connected with the 
modulated instability. 

Although the D3-D7' model also features a stable quantum Hall state, it was not 



included in the parameter space explored in 22,23 . Here we see for the first time 



how stability is recovered as the gapped, quantum Hall state is approached: rather 
than the range of unstable momenta shrinking, whole range runs off to infinity as the 
mass gap opens. 

We begin in Sec. [2] by reviewing the D2-D8' model, its construction and back- 
ground solutions. Then, in Sec. [3] we set up the fluctuation analysis of the metallic 
phase. Sec. |4]contains our results on the quasinormal mode spectrum and instabilities, 
and we conclude in Sec. [U 



2 Review of the D2-D8' model 



The D2-D8' system, introduced in 14 , consists of a probe D8-brane in the near- 
horizon background of iV near-extremal D2-branes, such that in the flat-space limit 
#ND = 6. We consider the limit of a large number of D2-branes N ^> 1 and large 
't Hooft coupling g s N ^> 1. Supersymmetry is completely broken, and internal flux 
is needed for stability. At weak coupling, the low-energy spectrum of bifundamental 
strings in a #ND = 6 system contains only charged fermions and no charged bosons; 
as a result, the holographic field theory dual is 2+1-dimensional SYM coupled to iV 
species of charged fermions. Similar constructions have been used in 2+1-dimensional 
D3-D7 intersections [2 16 . 
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2.1 Set-up 

The ten-dimensional near-horizon metric for a thermal D2-brane is 

(T 2 \ 
-^-y+M 2 rffi| 6 J , (1) 

where L 5 = 6ir 2 g s Nlg and the thermal factor f(u) = 1 — (^f) 5 - Here we use di- 
mensionless (lowercase) coordinates, which are related to the physical dimensionful 
(uppercase) coordinates as: x M = X^/L, u — U/L. In general, lowercase Latin letters 
will denote dimensionless quantities, and uppercase letters will denote the correspond- 
ing physical quantities. However, for the D8-brane worldvolume gauge field A^, we 
will define a dimensionless gauge field = A* ■ 

There is a black hole horizon at u = ut, and the corresponding Hawking temper- 
ature is T = ^-«y 2 . The metric on the S 6 is given by 

dQ 2 s6 = dip 2 + sin 2 ip (dB\ + sin 2 M0 2 ) + cos 2 ip (d£ 2 + sin 2 idQ\ + sin 2 f sin 2 9 2 d(p 2 2 ) , 

(2) 

where ip ranges from to tt/2, £, 9i, and 6> 2 range from to n, and 0i and 02 range 
from to 2n. The background dilaton and RR 5-form potential are 

e* - 9 ,(^f (3) 

C (5) = c(^)L 5 dn S 2 A dn s3 , (4) 

where c{jp) = | (sin-0 — | sin (3^) — ^ sin(5'0)). 

In this background, we consider a D8-brane probe filling the spacetime directions 
t, x, y, and the radial direction u. In addition, the D8-brane wraps an S* 2 x S 3 fibered 
over an interval in the internal S 6 . We will take the ansatz that the embedding in 
the ^-direction is a function only of u. The induced metric on the D8-brane is 

L- 2 ds 2 m = u^(-fdt 2 + dx 2 + dy 2 )+u-^(^ + u 2 ^du 2 

+-u~2 sin 2 ?f) dVt\ + -u _ 2 cos 2 ip dVt\ , (5) 



where ' = d u . 

We will also turn on the following three components of the D8-brane worldvolume 
gauge field: 

2na'F X y = h (6) 
2ira'F u0 = a' (u) (7) 
2-koIFq^ = L 2 6sin0! . (8) 

The constant spatial field strength corresponds to a background magnetic field h, 
the radial electric field is holographically dual to a charge density, and the magnetic 
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field wrapping the internal S 2 is needed to stabilize the embedding. Furthermore, we 
choose to work in the radial gauge, where a u = 0. 
The DBI action for the D8-brane is 

Sdbi = -/is J d 9 x - det(^ + Wi^) (9) 

= -N J duuWAGH , (10) 

where M = 8ir 3 T 8 v 3L 9 and v% is the dimensionless volume of spacetime. The functions 
A{u), G(u), and H{u) are defined as 

A(u) = l + /«V 2 -< 2 (11) 
G{u) = cos 6 ^ (6 2 M + sinV) (12) 

H(u) = 1 + ^ . (13) 



u 5 



The Chern-Simons action has a nonzero term due to the charge density and mag- 
netic field: 

S cs = -JL(2ira') 2 J C (5) AFAF = A/" y duc(tp)ha' . (14) 
2.2 Background solutions 



The D8-brane action (10) and (14) depends on clq only through its derivative. We 



can then integrate the equation of motion for ao to obtain 



jHd = d — hc(ip) = d(u) , (15) 



where we have defined 



'/(") /"' 2 \/^ ( 16 ) 



and where d is the constant of integration and d(u) is the radial displacement field. 
In terms of the boundary theory, d is the total charge density, while d(u) is the charge 
due to sources located in the bulk at radial positions below u. 
We can solve for a' in terms of d, obtaining 



On = d\ I - + ^ U — . (17) 
il d 2 + u 5 GH 



We then rewrite (16) as 



/ d 2 + u 5 GH 
HV 1 + fu 2 ^' 2 
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The ip equation of motion is given by 



2g 



d^G — ha' d^c(ip). 



(19) 



We can solve ( 19 ) at large u and find that ip — > in the UV. We identify the coefficient 



of the leading term in the asymptotic expansion of ip as the fermion mass and the 
coefficient of the sub leading term as the chiral condensate: 



m 
u 



ir 



+ 0(u 



(20) 



There are two classes of D8-brane embeddings with different IR behavior: black 
hole (BH) embeddings, which enter the horizon at ut, and Minkowski (MN) embed- 
dings, which smoothly cap off at some «o > «t as the wrapped S 3 shrinks to zero size. 
The DC conductivity can be computed straigthforwardly using the Karch-O'Bannon 
method [29] in the black hole phase and a generalization thereof [2] for the Minkowski 
phasej^] As argued in 14 , the MN embeddings holographically reproduce the prop- 
erties of a gapped, quantum Hall fluid. In our previous paper 19 , we analyzed the 
fluctuations of the MN solutions. 

Here we will focus on the quasi-normal mode excitations of the BH embeddings. 
The BH solutions correspond to an ungapped, metallic phase known as a holographic 
quantum liquid, with several unusual properties. At low temperature and no magnetic 

, T 5/3 



field, the longitudinal resistivity scales as p xx ~ T ' 14,31 , rather than like T as for 
a Fermi liquid. Such a fractional scaling exponent is due to hyperscaling violation; in 
a 2+1-dimensional system, thermodynamic quantities scale with an effective spatial 
dimension 2 — 6, where 9 is the hyperscaling violation exponent. For a Dp-brane 



background, 6 



(P-3) 2 



which for D2-branes gives 6 
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No concrete evidence 



of a Fermi surface has been seen for this, or any other Dp-Dq brane system (7 11 
though it has been argued that some sort of hidden fractionalized Fermi surface 
exists [6]. 

The MN embeddings exist only when d = meaning that the quantum Hall 
state is at filling fraction v 



1. In 14 



we found two types of transitions from the 
quantum Hall state to the metallic state. As the temperature is increased, there is a 
first-order phase transition from the MN embedding to a BH embedding. It persists 
for d/h Rj | as a first-order transition between two different BH embeddings. In this 
paper, we work at sufficiently low temperature such that we will completely avoid 
this transition. 

The second kind of transition occurs when the charge or magnetic field is changed 
so that d^\h. This causes the D8-brane to transition from a MN embedding to a BH 



embedding. We found in 14 that this process is not symmetric: increasing d leads to 



Alternately, one can use either the Kubo formula or the open string metric method 30 . 
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a continuous phase transition to a BH state J^] but decreasing d is discontinuous. There 
is a narrow region with d < |/i where the charge susceptibility < implying that 
the system is thermodynamically unstable]^] 

We will study further the nature of the transition from quantum Hall to metal- 
lic phases by examining the spectrum of fluctuations. In particular, will see this 



instability again in Sec. 4.1 as a diffusion mode with a negative diffusion constant. 



2.3 Temperature- independent scaling 

For the numerical analysis of the fluctuations, it is convenient to transform to a 
compact radial coordinate p defined as 

P= U -f- (21) 
u 

The boundary is at p = 0, and the horizon is now at p = 1. In this coordinate, the 
location of the horizon, and thus the temperature, has been scaled out. 

We can completely scale out the temperature by defining scaled quantities as 
follows 



m 


= m Urp 1 


(22) 


d 


_ 5 

= d Urp 2 


(23) 


h 


_ 5 

= k Urp 2 


(24) 


b 


1 

= b My 


(25) 


a 


= a Q Urr\ l 


(26) 
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_ 5 

= 9 Urp 2 


(27) 



and also rescaling the spacetime coordinates as (i,x,y) = u^ 2 (t,x,y). Using these 
scaled quantities, the action becomes 

S = -AfuJ/ 2 J dp (p- g / 2 VAGH + c(^)H) , (28) 



where now ' = d p and 



A(p) = l + j>y 2 -Ao 2 (29) 
G(p) = cos 6 i) (b 2 p + sin 4 A (30) 

H{p) = l + p 5 h 2 . (31) 



5 In the supersymmetric D3-D7 model, the careful analysis of 32 found that the analogous 
transition from MN to BH embeddings occurring at zero charge density was in fact third order. 
For the D2-D8' system, we have numerically verified that the continuous phase transition on the 
increasing d side of the MN state is also 3rd order. 

6 Viewing the systems in the grand canonical ensemble, as the chemical potential is lowered, there 
is a first-order phase transition where the d jumps discontinuously down from 
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The scaled equations of motion are 



a'o = -^Ljd (32) 



P 2 gH 
2gp 7 



d p (gHiP') = -L-d^G + ha'od^c . (33) 



3 Fluctuations 



3.1 Expansion 

We now allow for fluctuations of the following worldvolume fields: 

V> = $(p) + e6i/>(i,x,y,p) (34) 

ao = cLq(p) + e8a (i,x,y,p) (35) 

a x = e5a x (i,x,y,p) (36) 

d y = hx + eSd y (t,x,y, p) (37) 

a p = e5a p (i,x,y,p) , (38) 

with e a small parameter. We do not consider fluctuations in the internal spheres or 
dependence on the internal S 2 or S 3 coordinates, thus restricting to the s-wave. 
The Lagrangian can then be expanded to second order in e: 

£ = £ + e£i + e 2 £ 2 + ... • (39) 

Here Co corresponds to the background Lagrangian, and C\ vanishes on-shell by the 
background equations of motion. We derive the equations of motion for Sip, 5do, Sa x , 
and Sciy from £ 2 . In addition, since we are working in the radial gauge d p = 0, the 
equation of motion for 8d p becomes a constraint to be imposed on the fluctuations. 
We make the following wavelike ansatz for the fluctuations: 

6i/>(i,x,p) = 5i>{p)e-™ i+i% * (40) 
5a p (t,x,p) = 5d,(p)e- iU+ik£ , (41) 

where, because of the rotational symmetry in the (x, y)-plane, we can choose the 
momentum to be in the x-direction without loss of generality. 

We must be careful to consider only physical fluctuations rather than gauge arti- 
facts. To ensure this, we work with the fluctuation of the longitudinal electric field: 

Se x = k5a + u5a x (42) 

which is gauge invariant. As long as we start with initial conditions in radial gauge 
at some p which satisfy the constraint given by the equation of motion for 8a p , the 
constraint and gauge condition will be satisfied by the solution at all p. 
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Note that we have introduced a rescaled frequency and momentum which are 
related to their unsealed counterparts as (Co, k) = u T 3 ^ 2 (u, k). As a result, the longi- 
tudinal electric field then scales as 5e x = 5e x u T 5 ^ 2 . 

We finally obtain a coupled system of linear, second order equations for the fluc- 
tuations. The 5tp equation of motion is 



d P ( X \^f + (i + A?) <^') J + w P ^ GH / d ^ G f^ 



-i5a y d p [ xj^fhki)' 



ikba y I djca - W — h I - dip I a d^c 

- k 2 (l - p 4 Pa' 2 ) 



Gp 3 (HCo 2 u2( , M ^ l2 \\ I AGE ( d$ log G\ 2 ^ /Atf^G 



AH \ f \ ' U J J Vp 9 V 2 J V Gp 9 2 



■jjjkt/j a 5e x - 5a y \hd^c - J — a ) . (Li) 

The equation of motion for 5a is 



^ ( c + V % h <) - x = \l% ta '^ + - f Xl%( l + fp2 ^ Si * ' (44) 

where we have defined the function J as 



A 3 

The equation for 5a x is 



• (45) 



c+ W— /ia„ + <9 P 



A# u / My AHp cu 



^ka' o m + (1 + /p 2 ^' 2 ) 5e~, (46) 



and for Sa y is 



A (-< + J^^h) H + ^ - (1 + fm fj «, (47) 
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Finally, the equation derived by varying Sa p is a constraint which imposes the gauge 
condition a p = 0: 




, TT —Sa' 
AHpuo x 



(4£ 



The constraint can be used to eliminate one of the fluctuations in terms of the others. 



Note that using the constraint, the 6qq equation (44) and the Sa x equation (46) are 



identical, implying that there are only three independent equations of motion. 



3.2 Numerical methodology 

Our goal is to solve the equations of motion for the fluctuations and find the normal- 
izable modes. The condition for normalizability is that the fluctuations vanish at the 
boundary: 

5e x (0) = 5a y (0) = <ty(0) = . (49) 

At the horizon, the fluctuations have infalling boundary conditions. 

In practice, to numerically solve the coupled equations for the fluctuations at 
fixed values of Cj and k, we vary the initial conditions at the tip and shoot out to the 
boundary to look for solutions which are normalizable. We must solve three coupled, 
linear, second-order differential equations for three independent functions. 



An efficient technique for finding such solutions is the determinant method 21 33 
Rather than attempting to scan, for each Cj and k, the full three-dimensional space 
of initial conditions, we choose a basis of linearly-independent initial conditions, 

(8p(l), 5e x (l), Sa y (l)) = {(1, 0, 0), (0, 1, 1), (0, 1, -1)} , (50) 

for which we numerically solve the equations of motion. We then compute the fol- 
lowing determinant of boundary values: 



Sip 1 54> n 54> HI 
(let I 5e{ del 1 b~e{ u 
b~a\ 5^ b~a\ u 



(51) 



where the superscripts refer to the basis elements of (50). If the determinant vanishes, 
there is at least one linear combination of the basis vectors for which 5ip(0), Se x (0), 
and Sa y (0) simultaneously vanish, signaling the existence of a normalizable mode. We 
repeat this procedure for various values of Cj and k; for k = 0, a scan in Cj yields the 
spectrum, and allowing k > enables us to compute the dispersions. 

For the figures presented in Sec. |4| we have chosen particular values of the pa- 



rameters ut, m, d, h, and b for illustrative purposes. As discussed in Sec. 2.3, the 



system is invariant under a scaling of the coordinates and parameters so there is, in 
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fact, only a four-dimensional parameter space of BH solutions. We used this rescaling 
symmetry to eliminate ut from the equations. However, it makes more physical sense 
to work at fixed mass m and background flux 6, and vary the exogenous variables ut, 
h, and d. With the exception of Fig[l3j which shows the dependence of the instability 
on the mass, all of the plots in Sec. [4] are at fixed m = 1 and 6=1. 



4 Numerical results 

4.1 Hydrodynamical diffusion mode 

The dominant quasinormal mode is the excitation with the smallest imaginary part, 
which then decays the slowest. At long wavelengths, that is k — > 0, we find the longest 
lived mode is a purely imaginary hydrodynamical diffusion mode, with dispersion 



-iAk 2 + 0(k 3 



(52) 



as shown in Fig. [T} This type of diffusion mode is generically found in systems with 
charge conservation; a local charge excess can only dissipate by diffusing through the 
medium. Such a mode has been seen in other similar brane models, such as D3-D7 



and D3-D7' 21,22,24 



Im cj 




Figure 1: The dispersion of the hydrodynamical mode at h = (left) and h = 0.35 
(right) for ut = 0.3, d = 0.2, m = 1, and 6 = 1. The curve is a fit to the quadratic 
dispersion ([521); from this fit, for h = 0, A = 3.83, and for h = 0.35, A = -4.96 x 10~ 3 . 



A negative A indicates that the system is unstable to charge density fluctuations. 



When the magnetic field is turned off, the diffusion constant grows with charge, 
as shown in Fig. [2] (left). In a nonzero magnetic field, Fig. [2] (right), A similarly grows 
for d > At exactly d = |/t, A = 0; the system is in the gapped quantum Hall 
state, and the diffusion mode disappears. 

At vanishing charge density, as the magnetic field is increased, A initially increases 
to a maximum and then decreases asymptotically to zero, as illustrated in Fig. [3] (left). 
At nonzero d, for some h sufficiently larger than |d, A similarly reaches a maximum 
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io : 



0.15 - 



Figure 2: The diffusion constant A of the hydrodynamical mode as a function of 
charge density d, (left) at zero magnetic field h and (right) at magnetic field h = 0.15, 
with ut = 0.3, m=l, and 6=1. The solid blue portion of the curve indicates that 
the system is stable. The black section where A < is unstable. The dashed red 



regions are subject to a modulated instability; see Sec. 4.3 



and then decays to zero; see Fig. [3] (right), although the decay is not shown. For h 
below |d, A is decreasing. 

However, the most striking features of both Figs. [2] and [3] are the regions, drawn in 
black, where A < 0, as illustrated in Fig. [T] (right); this implies that Im uj > and the 
diffusion mode is unstable. Physically, this instability means that long wavelength 
fluctuations in the local density grow rather than diffuse away. This is actually just 



another indication of the thermodynamic instability discussed in Sec. 2.2 In fact, as 
can be seen in Fig. |1J the narrow region adjacent to the quantum Hall state where 
A < matches precisely with the region where |nf- = |^ < 0. 



4.2 Zero sound 

A generic consequence of a Fermi surface is the existence of a zero sound mode, 
which is a collective excitation of the Fermi surface. Landau-Fermi liquid theory 
makes precise predictions for the properties of the zero sound. At zero temperature, 
the dispersion is 



uj 



(53) 



i , where d is the number of spatial dimensions. The 

2 



The zero sound speed c s = 

quantum attenuation T q ~ k 2 /fi is due to decay into quasiparticle-quasihole pairs 
At nonzero temperature, thermal fluctuations further attenuate the zero mode. For 



kfj, ^> T the system is collisionless, and the zero sound with dispersion (53) persists 



although T q is temperature-dependent. When k/j C T 2 , the system is hydrodynami- 
cal, and the zero mode is completely damped and disappears. 

Despite the lack of direct evidence for a Fermi surface, a holographic zero sound, 
which looks very much like the usual Fermi liquid zero sound, has been identified in a 
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Figure 3: The diffusion constant A plotted versus magnetic field h, (left) at d = 
and (right) at d — 0.2, with ut = 0.3, m — 1, and 6=1. The solid blue part is stable, 
the black part is unstable, and the dashed red part is unstable, but only at nonzero 



momentum; see Sec. 4.3 



variety of Dp-Dq systems |8}|TT||20||22}]24||34| . The D2-D8' model is no exception, and 
the zero sound here shares many similarities with those of other holographic models. 
In Fig. |5j we plot the imaginary part of the dispersion for the longest lived modes, at 
zero magnetic field. Fig. [6] shows the dispersions of these low lying modes for various 
values of the magnetic field. 

At nonzero temperature and for small or vanishing magnetic field, the zero sound 
first appears at the transition from the hydrodynamical regime to the collisionless 
regime. At that point, the purely dissipative hydrodynamical mode merges with 
another purely imaginary mode, see Fig. [6] (left). The two modes then go off the 
imaginary axis into the complex plane to form the two branches of the zero sound 
with a nearly linear dispersion. The location of the crossover from the hydrodynamical 
to collisionless regime is temperature-dependent: as the temperature is reduced, it 
occurs at smaller k, approaching k = as the temperature goes to zero. 

At small magnetic field, this second imaginary mode merges at a small k with 
a third purely imaginary mode, forming another pair of complex modes. For these 
long-wavelength modes, u is nearly independent of k and scales linearly with magnetic 
field. 

As the magnetic field is increased, this second merger occurs at even larger k, 
see Fig. [6] (middle). Eventually, at a critical magnetic field h m , the two mergers 
meet, leaving a single, purely imaginary diffusion mode and a pair of complex modes 
(Fig. [6] (right)); we can identify the complex modes as a massive zero sound, whose 
mass grows linearly with h, at least for small h, as shown in Fig. [7} 
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Figure 4: The unstable regions in the (/i,d)-plane with ut = 0.3, m = 1, and 6=1. 
(Left) Overall plot showing all phases. (Right) A zoom into the region of small d and 
h. The straight black line, h — Id is the stable QH state. The region above that, best 
seen in the right hand plot and labeled "Unstable", is the region of thermodynamic 
instability. The red curve is where k m i n hits zero. On the green curve A crosses zero, 
and at the blue curves k min = k max . The dashed black curve with open boxes denotes 
the boundary of the region where |^ < and coincides to a very high accuracy with 
the red and green curves. 



4.3 Modulated instabilities 

Following the lowest purely imaginary mode to larger k, in some cases, Imu becomes 
positive for a range of momenta, signaling the onset of an instability; a typical such 
dispersion is shown in Fig. |8j Because the unstable modes are at k > 0, this decay 
leads to spatial modulation, breaking translation invariance. 



This modulated instability is a variant of the type first investigated in 26 and 



subsequently seen in many other holographic contexts 22 23 27 . Chern-Simons gauge 
theories are prone to modulated instabilities in the presence of electric fields but are 
typically stabilized by magnetic fields. In many holographic contexts, the radial 
electric field dual to the charge density is responsible for the instability. Here, the 
situation is somewhat more complicated, as the fluctuations are completely coupled 
at k > 0, and the instability can persist down to zero density. 

Roughly speaking, we find the instability occurs at low temperatures, large mass, 
and in the vicinity of the quantum Hall state. That is, for fixed h and d, there is 
a critical temperature u c T above which the instability disappears, and sufficiently far 
from d = the system is stable. In addition, the unstable region of parameter 
space shrinks as the mass is reduced. 

We can illustrate this for representative parameters in various ways. Fig. [4] shows 
the region of the (h, d)-plane, at fixed temperature, where the modulated instability 



occurs. Fig. [9] shows function of d for two different values of h, while Fig. RL0 

plots function of h for two different values of d. 
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Figure 5: The imaginary parts of the dispersions for the longest lived quasi-normal 
modes, for d = 0.2, h — 0, ut — 0.3, m = 1, and 6=1. Purely imaginary modes 
are denoted by solid lines, while dashed curves represent complex modes. The region 
near the origin is shown in more detail in Fig. [6] (left). The hydrodynamical diffusion 
mode, seen in Fig. [I] (left), is the solid black curve. The dashed red mode is the zero 
sound. The solid magenta mode crosses the fc-axis at k — 0.91, signaling the onset of 



a spatially modulated instability; see Sec. 4.3 



rp IS decreasing as h increases, and, as expected, the 

2. 



The behavior of the critical temperature is not monotonic either as a function of 
density or magnetic field, but, surprisingly, it depends on which side of the quantum 
Hall state one is on. For h > |d, u 

magnetic field is stabilizing the system. However, for h < |d, increasing h raises Uj<, 
evidently destabilizing the system. Note, in particular, that for d > increasing 
the density decreases the critical temperature. This is in contrast to what was seen 
in the D3-D7' model 22 , where goes down with increasing d. 

We present in Figs. 11 and 12 the range in k of unstable modes for four slices 
across the (h, <i)-plane shown in Fig. |4j Stability away from d = §/t occurs when 
kmin = k max . The k m in {k ma x) is the minimum (maximum) momentum above (below) 
of which there is an unstable mode. Both Figs. 11 and 12 include domains where 
kmin = 0; these correspond to regions of homogeneous thermodynamic instability. 

Although there are unstable regions on either side of d = the quantum Hall 



state itself is stable 19 



10 



Counterintuitively, as can be seen from Figs. |9| and Fig 
the critical temperature increases as the quantum Hall state is approached, implying 
that the system is becoming more unstable, rather than less. However, as is evi- 
dent from Figs 



11 



and 



12 



\h is 



both k min and k max grow toward infinity as d 
approached. The exception is when the quantum Hall state is approached from the 
thermodynamically unstable region, in which case k m i n stays zero. 

Finally, we illustrate the effect of the mass on stability in Fig. 13 we again plot 
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Figure 6: The dispersions of the lowest lying quasi- normal modes at h = (left), 
0.005 (middle), and 0.01 (right), with d = 0.2, u T = 0.3, m = 1, and b = 1. The 
real part of u; is shown by dashed curves, and the imaginary part by solid curves. 
The purely dissipative hydrodynamic mode is black, and the zero sound is red. (Left 
and middle) For h < h m « 0.01, the hydrodynamic mode merges with a second, 
purely imaginary blue mode and then goes into the complex plane to become the 
zero sound. The blue mode merges at lower k with a magenta purely imaginary 
mode, becoming two orange, long-wavelength, complex modes with nearly constant 
dispersion. (Middle) This second merger point moves to larger k as h is increased. 
(Right) Above the critical magnetic field h > h m , the two merger points have met, 
leaving only the dissipative mode and the now massive zero sound. 



the range in k of unstable modes, this time as a function of m. The range of unstable 
momenta shrinks as the mass is reduced. For a given d, h, and ut, there is a sufficiently 
small m such that the modulated instability disappears entirely. 



5 Discussion 

Our investigation of the fluctuations of the D2-D8' model found a quasi-normal mode 
spectrum in many ways qualitatively in line with those of other related probe brane 
systems. We found a long- wavelength diffusion mode in the hydrodynamical regime 
and a holographic zero sound in the collisionless regime. An applied magnetic field 
modifies the spectrum and, in particular, generates a mass for the zero sound. At 
low temperature and for the charge density and magnetic field in a range near the 
quantum Hall state, the homogeneous phase suffers a spatially modulated instability. 
In certain parameter regions the minimal spatial modulation goes to zero continuously, 



reproducing the thermodynamic instability found in 14 . 

We did, however, find some significant differences, in particular, with D3-D7' 
model to which the D2-D8' model is most closely related. Most notably, we found 
here an unstable region near the quantum Hall state, while in the D3-D7' model, the 
instability was found only at high density and low magnetic fields. In addition, we 
found smaller masses led to improved stability, which is exactly the opposite of what 
was seen in the D3-D7' case. Considering that the instabilities stem from essentially 
the same mechanism, this is quite surprising. One potentially related issue is that 
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Figure 7: The real (dashed red) and imaginary (solid blue) parts of u for the zero 
sound mode at k = as a function of h, with d = 0.2, ut = 0.3, m = 1, and 6 = 1. 
The mass, given by the real part, increases linearly with magnetic field for h < 0.3. 
Notice that the region, where the merging of the modes depicted in Fig. [6] occurs, lies 
very close to the origin in this plot. 



in 



22 23 , the two independent fluxes of the D3-D7' model were chosen to be equal, 
which excluded the possibility of a quantum Hall state. If one of the fluxes were shut 
off, the unstable region of the D3-D7' could possibly also be centered around the now 
allowed quantum Hall state. 

An important open problem is to find the ground state to which this modulated 
instability ultimately leads. As the unstable modes have nonzero k, it is reasonable 
to expect that translation invariance will be broken in the ground state. Modulated 
states, such as spin and charge density waves, with spontaneously broken translation 
symmetry, appear in various condensed matter contexts^] Explicit modulated ground 
states have been studied in some simpler holographic systems exhibiting modulated 
instabilities |27|. In this model, finding such a state is technically more challenging 
and is a subject of ongoing investigation. Because all the fields are coupled, we 
anticipate a striped phase in which all quantities, such as charge density, currents, 
and the fermion condensate, are spatially modulated. 

As the unstable region was found in the vicinity of the QH state in this model, 
it would be interesting to check if the quasi-normal modes can be matched to the 
normal modes of the QH state in this region. In particular, one could expect that 
the zero sound approaches the lowest QH mode, which is a magneto-roton in a vast 
region of the parameter phase 19 28 , as the magnetic field approaches the locking 
value of the QH phase. As the numerical analysis of the quasi-normal modes close to 
this phase is very involved, our attempts to clarify this issue remain inconclusive. 



7 For reviews of spin and charge density waves in condensed matter systems, see 
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Figure 8: The dispersion of the lowest purely imaginary mode with d = 0, h = 0.1, 
ut = 0.3, m — 1, and 6=1. Between k min = 0.712 and k max = 2.189, Imw > 
signaling the existence of a modulated instability. 
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Figure 9: The critical temperature ut below which the system suffers a modulated 
instability as a function of d, (left) for h = and (right) for h = 0.3, with m — 1, 
and 6=1. 



Another outstanding issue is the existence of a Fermi surface. Fermi surfaces have 
been argued to exist in holographic quantum liquids, but no compelling evidence 
for them in probe brane models has been found. It may be that if the charged 
fermions are holographically encoded in bulk degrees of freedom at the black hole 
horizon, the Fermi surface is fractionalized and difficult to detect by the usual gauge 
invariant probes. For this reason, it may be useful to focus on those models with 
charged fermions dual to sources outside horizon. In brane models with #ND = 6, 
background magnetic fields can induce, via Chern-Simons terms, charge along the 



holographic radial direction. This effect has been studied in the D2-D8' 14 , D3 



D7' |2|, and Sakai-Sugimoto models 36 , where the charge carried on the horizon 



and the charge spread radially along the brane behave as two fluids with distinct 
properties. Perhaps, these can serve roughly as probe brane analogs of the electron 
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star 37,38 



We end this paper by discussing another possible outgrowth of this work. In the 



D3-D7' model, in 39 , several non-trivial examples has been uncovered where the 
weak-coupling results for certain observables has been matched with the results of 
holographic strong-coupling calculations. It would be interesting to study whether 
this analysis could be carried over also for the D2-D8' model. However, this would be 
much more involved in the present context due to background being only conformal 
to AdS and the lack of massless probe embeddings. 
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Figure 11: The minimum (solid blue) and maximum (dashed red) momentum k where 
the lowest purely imaginary mode is unstable as a function of charge density d, (left) 
at h = and (right) at h = 0.15, with ut = 0.3, m = 1, and 6 = 1. Note that for 
h = 0.15 (right), k min = for d G (0.065,0.1), implying that there is a homogeneous 
instability. In the right hand plot, because we have zoomed in to show the behavior 
around d ~ 0.1, the maximum momentum is above the region shown; for d G (0.1, 0.2), 
> 25 

max r~~j 
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